Directed current in the Holstein system 
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We propose a mechanism to rectify charge transport in the semiclassical Holstein model. It is 
shown that localised initial conditions, associated with a polaron solution, in conjunction with a 
nonreversion symmetric static electron on-site potential constitute minimal prerequisites for the 
emergence of a directed current in the underlying periodic lattice system. In particular, we demon- 
strate that for unbiased spatially localised initial conditions, violation of parity prevents the exis- 
tence of pairs of counter-propagating trajectories, thus allowing for a directed current despite the 
time-reversibility of the equations of motion. Occurrence of long-range coherent charge transport is 
demonstrated. 

PACS numbers: 05.60.Cd, 05.45.Ac, 05.60.-k, 05.45.Pq 



Since the pioneering work of Landau and Pekar it is 
established that when an electron in a lattice interacts 
locally with the phonons it can become self-trapped by 
the lattice distortion it creates [H,0]- The generated 
quasi-particle consisting of the localised electron and its 
associated local lattice deformation is called a polaron 
(an electron surrounded by a phonon cloud). Frohlich 
(large polarons) and Holstein (small polarons) elaborated 
significantly on the consequences of the polaron concept 
[3| , 0] • Later on this idea of a localised quasi-particle was 
used by Davydov to propose a model for coherent energy 
and/or charge transport in macromolecules 

Recently the emergence of a directed current trig- 
gered by an external time-dependent field with zero 
mean which can be of stochastic or deterministic na- 
ture has attracted considerable interest. Applications 
include charge transport in semiconductor heterostruc- 
tures 1101 , current generation in semiconductor superlat- 
tices [ll[, nano-engines [13], directed collective energy 
current in spatially extended systems [l3j], superconduc- 
tors 14 1, spin transport [Hj], and the motion of cold 
atoms in driven optical potentials [ill ]. 

The necessary conditions for rectification of the cur- 
rent, based on symmetry investigations of the external 
field and the underlyin g st atic potential, have been pre- 
sented in 11 j , [I?) , and [lH . To be precise, all symmetries 
that, to each trajectory, generate a counterpart moving 
in the opposite direction, need to be broken. Further- 
more, the phase space has to be mixed, with coexisting 
regular and chaotic dynamics [l9j . This is achievable by 
imposing a time-dependent external force that is peri- 
odic but not symmetric under time reversal [l3] — El| - I n 
extended chaotic systems a nonzero current can be ob- 
tained as the time-averaged velocity of an ensemble of 
trajectories in the chaotic component of phase space and 
the chaotic transport proceeds ballistically and direct- 
edly [llj], [2l|. Recently the emergence of directed flow 
in autonomous Hamiltonian systems has been demon- 
strated in [24]. There it is shown that internal parity 
violation and localised initial conditions, which are unbi- 
ased, in conjunction with transient chaos conspire to form 
the physical mechanism for the occurrence of a current. 

In this manuscript we propose a symmetry-breaking 



mechanism to accomplish directed electron (charge) 
transport in the context of the semiclassical Holstein sys- 
tem. The Hamiltonian of the semiclassical Holstein sys- 
tem consists of three parts 
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The first term 



E 



{U n \ Cn \ 2 - V ( C * n C n ^ + CnC^)} . (2) 



describes quantum-mechanically the motion of an ex- 
cess electron in the tight-binding approximation over the 
units of a one-dimensional chain where c„ determines 
the probability amplitude to find the electron residing 
at site n. V is the transfer matrix element (its value is 
determined by an overlap integral) being responsible for 
the nearest-neighbour transport of the electron along the 
chain and U n denotes the on-site potential (on-site en- 
ergy of the electron). The second term represents 
the Hamiltonian for the classical dynamics of (local) Ein- 
stein oscillators 
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(3) 



The coordinates q n quantify the displacements from equi- 
librium position of the oscillator at site n, p n is the cor- 
responding canonically conjugate momentum, loq is the 
harmonic vibrational frequency, and m is the mass of an 
oscillator. The interaction between the electronic and the 
vibrational degrees of freedom is due to the modifications 
of the electronic on-site potential U n by displacements of 
the oscillators from their equilibrium positions and is is 
described by the following interaction Hamiltonian 



Hj, 
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where a regulates the electron-oscillator coupling 
strength. 

For a dimensionless representation we introduce the 
following time scale: t = ljq t. The dimensionless repre- 
sentation of the remaining variables and the parameters 



2 



of the system follows from the relations: 
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In what follows we drop the tildes. 

The equations of motion derived from the Hamiltonian 
given in Eqs. (J2|)-(j4|) read as 

dc n 

ir—^- = U n c n - (c n+ i + Cn-i) + aq n c n (6) 
at 



d 2 g n 
dt 2 



a\c n Y 



(7) 



The parameter r = Hujo/V, appearing in the l.h.s. of 
Eq. ©, determines the degree of time-scale separation 
between the (slow) electronic and (fast) vibrational pro- 
cesses. 

In the following we construct standing localised solu- 
tions. To this end we consider the corresponding sta- 
tionary system. For the oscillator system the condition 
ijn = yields the stationary displacements 
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(8) 



from which, after substitution in Eq. ^ , one obtains the 
following discrete nonlinear Schrodinger equation 



■ ^IL - TT 
dt 



(C n +l + Cn-l) - "f\c n \ 2 Cr, 



(9) 



with 7 = a 2 . Substituting the ansatz c n (t) — 
4> n exp(—iE/Tt) with <f> n 6 R in Eq. © one arrives at 
the system of coupled difference equations 



E<t>n = U n 4> n 



i+l + 4>n-l) 



(10) 



This system supports localised solutions corresponding 
to standing localised electrons. Following the methods 
outlined in Refs. [22j],[23| we derive the profile of the 
standing localised electron wave function. Subsequently, 
with the help of Eq. (jHJ one obtains the corresponding 
localised lattice oscillators forming in conjunction with 
the localised electron the polaron compound. One utilises 
that polaron solutions are obtained as the attractors of 
the map 



{4>}^y>}=H{cj>}/\\H{ct>}\\ 



(11) 



where the operator H is determined by the right hand 
side of Eq. (fT0|) and the norm of the state H{cf>} is de- 
fined as \\H{(j)}\\ = VEnO^W) 2 - The ma P itera- 
tions are started with a completely localised state, i.e. 
{</>n°^} = <5 n ,n , and act on it with the operator H. After 
each application of H the resulting vector is normalised 
and the iteration procedure is terminated when conver- 
gence is attained yielding the polaron state of lowest en- 
ergy. In Fig. [1] we show the profile of the electron occu- 
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FIG. 1: Profile of the stationary electron occupation proba- 
bility displayed for a lattice segment around the central site 
no = 0. The parameter values are given by a = 1.25 and 
ljq = 2. Solid line: The symmetrical on-site potential with 
a = Uo = 0, Dashed line: The asymmetrical on-site potential 
with a = and U = 0.035. 



pation probability 0^ for two different situations. The 
first case is for uniform electron on-site potential U n = 
(solid line in Fig. [1]). (For U n = U = constant, a sim- 
ple gauge transformation c n (t) — c n (t) exp(iUt) removes 
the on-site term from the equation of motion.) The peak 
of the profile is located at the central site no = and 
the occupation probability decays to either side in an 
exponential fashion. Note the inversion symmetry with 
respect to this site, i.e. <j> n = <p- n . For the second case 
we chose a sawtooth-like profile of the electron on-site 
potential of the following form 



' Uo, n = 3k; 
r n -a=<2[/ , n = 3k+l; 

,3C/ , n = 3fc + 2. 



(12) 



where k = 0, ±1, ±2, .... Notice that an electronic on-site 
potential of the form in (fT2j) breaks the reflection sym- 
metry no + n no — n of the lattice system inducing 
in this way an asymmetric potential biasing motion to 
the right. The lattice system is still periodic exhibiting 
translational invariance. In fact this potential possesses 
the shortest possible 'unit cell', namely of period-3, vio- 
lating the reflection symmetry of the lattice system. The 
resulting asymmetry is reflected in the profile of the static 
polaron solution which for an asymmetric potential with 
a = and Uq — 0.035 is depicted by the dashed line in 
Fig. [TJ Remarkably the peak of the occupation probabil- 
ity is not situated at the central site no = but is shifted 
one lattice site to the left. As in the symmetrical case 
the occupation probability is exponentially localised at a 
single site. 

For later use we introduce here for a lattice comprising 



3 



N sites the participation number, P, denned as 



P 



1 



N 



(13) 



quantifying the spatial extension of the electron occu- 
pation probability (degree of electron localisation). An 
electron completely localised at a single lattice site corre- 
sponds to P = 1 whereas P — N corresponds to a com- 
pletely delocalised electron for a finite lattice of N sites. 
For the occupation probabilities displayed in Fig. [1] one 
obtains P ~ 7 for both the symmetrical and the asym- 
metrical case. 

For a quantitative assessment of directed transport we 
consider the current density J being proportional to 



J(t) 



1 N 



i 



N 



i-^Im[<(t) (c n+1 (t) 

71=1 



Cn-l(t))] • (14) 



In the following we discuss the symmetry features of the 
lattice system. We are in particular interested in the 
emergence of a directed current facilitated by the asym- 
metric potential. 

In order to obtain a non-zero current (net flow) , all spa- 
tiotemporal symmetry operations that leave the system 
of equations of motion invariant and invert the sign of the 
current density must be broken. Otherwise the current 
density J completely vanishes. Regarding the symme- 
tries of the semi-classical equations of motion we note 
that Eqs. ([6 ]) .([7 | are invariant under the time-reversal 
operation t i-> —t. For the canonically conjugate mo- 
mentum and coordinate variables of the classical system 
of Einstein oscillators this is equivalent to the trans- 
formation {{p n {t)},{q n {t)}) i-> ({-Pn(-t)}, {q n (-t)}). 
For the quantum discrete Schrodinger equation the time- 
reversal operation, t i— > —t, is connected with complex 
conjugation of the probability amplitudes, c„, which also 
reverses the sign of the current density. Passing from 
the pairs of complex- valued canonically conjugate vari- 
ables, {c ra ,c*}, to real- valued canonically conjugate mo- 
mentum and coordinate variables via the canonical trans- 
formation c„ = (Q n + iP n )/y2 and using the repre- 
sentation c„ = r n exp(i#„) one obtains the relations 
P n = y/2r n sm(i8 n ) and Q n = \[2r n cos(i9 n ). Upon 
time-reversal operation, t <— > —t, the corresponding inver- 
sion of the sign of the phases after complex conjugation, 
i.e. c* = r„ exp(— i6 n ), is equivalent to the transforma- 

tion ({P„(*)},{Qn(*)}) -> ({-P»(-t)},{Qn(-*)» ■ 

In general, solutions of the system are of the form 
X(t) = ({Pn(t)}, {q n (t)}, {<(<)}, {c n (t)}) ■ (15) 
Thus applying the time-reversal operator yields 

H{Pn(t)}Aqn(t)},{c;M},{Cn(t)}) 
= {(-Pn(-t)}, {?„(-*)}, {Cn(-t)}, {<(-*)}), (16) 



and hence, if X is a solution, then so is fX. As for the 
implication of time-reversibility symmetry with regard 
to the net flow, let a solution, starting from some initial 
condition A(0), be evolved in time up to a finite obser- 
vation time T at which the forward trajectory arrives at 
the point X(T) on the constant energy surface. Subse- 
quently, under application of the time-reversal operation, 
for the classical subsystem the signs of the momenta {p n } 
at this point, X(T), are reversed while the complex con- 
jugation is applied to the probability amplitudes {c„}. 
Letting then the solution evolve once again, with tX(T) 
as the initial condition, the corresponding backward tra- 
jectory traces back the path of the forward trajectory in 
coordinate space {{q n },{Qn} = Re({c„})). For a mi- 
crocanonical ensemble, the initial conditions X(0) and 
tX(T) are equally selected points from the constant en- 
ergy surface. Thus, for systems with time-reversibility 
symmetry and uniformly distributed initial conditions 
populating the whole energy surface there is no preferred 
direction of the flow thus preventing the emergence of a 
current. 

However, considering a static polaron solution as the 
initial state for the dynamics corresponds to initial condi- 
tions being exponentially localised at a specific site, no, 
of the lattice. Therefore, only a confined region of the 
energy surface contains initial conditions. Suppose the 
directed motion of a polaron can be instigated; then the 
polaron undergoing directed motion leaves the domain of 
localised initial conditions on the lattice (the extension of 
which is determined by the participation number, P, de- 
fined above) and contributes to a net current. Crucially, 
at the end of the observation time, T, the polaron (as- 
suming it travels directedly far enough while retaining its 
localised shape) lies outside the domain of the localised 
initial conditions. Consequently, the initial condition of 
the associated backward moving polaron, compensating 
the contribution of the forward moving polaron to the 
net current, is not contained in the set of localised initial 
conditions. We underline that this alone does not imply 
the emergence of a current in the system for such sets of 
initial conditions [24j |. 

In fact, in order to obtain a non-zero current density 
the remaining invariances of the system, which induce 
symmetries ensuring that there exist trajectories which 
mutually compensate each others contribution to the net 
flow, must be broken. 

In this context we note that the system exhibits reflec- 
tion symmetry n$ + n i— > no — n if and only if the elec- 
tronic on-site potential U n possesses reflection symmetry. 
In more detail, the following parity symmetry operations 
leave the system of equations of motion invariant and 
reverse at the same time the sign of J: 

n Q + n ^ n - n , U no+n (17) 

^no+n 1 ^ Cno—n 3 Qno+n 1 ^ Qno— n • (1^) 

As a consequence, in the presence of parity (reflection) 
symmetry there exist pairs of counterpropagating trajec- 
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tories, 

x +(t) = ({Pn +n(t)},U n +n Cn +n 

(*)}) 

(19) 

and 

*-(*) = ({-P»o-n(*)},{9no-n(*)},{< -„(t)},{c„ - n (*)}) 

(20) 

annihilating each others contribution to the net current 
density. 

Crucially, an asymmetrical electron potential such as 
the one of the form in (|12p breaks the parity symmetry 
and thus allows for the occurrence of a non-zero current 
density. 

So far we dealt with standing polaron solutions having 
zero associated kinetic energy. As a first step we initiate 
the motion of the polaron through suitable initial con- 
ditions of the momenta of the lattice oscillators Jp„(0)} 
targeted in the direction of the pinning mode |25j . To 
be precise we use the following initial conditions for the 
numerical integration of the system; 

{g n (0),p„(0),Re(c„(0)),Im(c n (0))} = {<?„, A£, </>„, 0} , 

. (21) 

with the normalised momentum part £ of the pinning 
mode and A regulating its amplitude. Here, we exploit 
the lowest frequency pinning mode although others of 
higher frequency could also be applied. We solve numer- 
ically the set of coupled equations using a Runge-Kutta 
method for a lattice consisting of N = 105 sites. The 
norm conservation ^ n |c n (i)| 2 = 1 (as well as the con- 
servation of energy) was monitored during the integration 
procedure to guarantee accurate computations. We em- 
ploy periodic boundary conditions. We stress the relation 
between the extension of the polaron and the length of 
the lattice P < N. 

The effect of asymmetry becomes apparent in Fig. [2] 
which displays the spatiotcmporal evolution of the elec- 
tron occ upat ion probability for pinning mode amplitudes 
A = ±V0.15. We underline that interchanging the sign 
of the amplitude A <H> — A entails merely the reflection 
anti-symmetry p„(0) <H- —p n (0) in the initial momenta 
under preservation of the kinetic energy, thus, imposing 
no bias. For the pattern shown in the top panel (positive 
A) the electron tends to move to the right in an initial 
phase. However, it soon becomes trapped at a lattice site. 
In stark contrast with this, its counterpart in the bottom 
panel (negative A) exhibits coherent motion to the right 
under maintenance of the localised pattern (despite the 
initial slight broadening of the profile). Strikingly, radi- 
ation losses due to the emission of phonons seem to be 
suppressed. In this case the electron is accordingly ac- 
companied by the associated localised pattern of lattice 
oscillators (not shown). 

Furthermore, we utilised ensembles of unbiased initial 
conditions for the momenta to instigate the motion of po- 
larons. To be precise, kinetic energy is provided locally to 




FIG. 2: Spatiotemporal evolution of the electron occupation 
probability (o.p.) for a lattice comprising N = 105 sites and 
periodic boundary conditions are imposed. The parameter 
values are given by a = 1.25, ujq = 2, r = 25, a = 0, and 
U = 0.035. Upper (lower) panel A > (A < 0). 

the system as follows: In the region of the lattice support- 
ing the localised standing-electron lattice-oscillator com- 
pound whose spatial extension is determined by the par- 
ticipation number P, kinetic energy of magnitude E^in 
is injected in such a way that the initial values of the 
momenta, {p„(0)}, are distributed uniformly on an iso- 
energetic sphere in momentum-space such that the rela- 
tion 

1 [P/2] 

Ekin = 2 2 p «o+«- ( 22 ) 

n=-[P/2] 

is fulfilled and [ • ] denotes the integer part of P/2. Notice 
the symmetry pi O —pi and i — uq — [P] , . . . , no + [P] . 
Hence there is no bias contained in the ensemble of initial 
conditions {p n (0)}. We emphasise however the localised 
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FIG. 3: Temporal behaviour of the mean current density 
(J(t)) for an ensemble of 1000 initial conditions distributed 
according to Eq. (122 p with Ekin — 0.1. The parameter values 
are given by a = 1.25, ojq = 2, r = 25, a = 0, and (To = 0.035. 

character of the initial conditions. For a quantitative 
assessment of the efficiency of the transport-rectifying 
mechanism we consider the mean current density (J) 
where (•) denotes the ensemble average. In Fig. [3] we 
show the temporal behaviour of the mean current density 
for our simulation for which the motion of the standing 
asymmetrical polaron shown in Fig. [T] is instigated with 
the inclusion of kinetic energy in the initial conditions in 
the lattice segment n £ [—4,4]. To this end an ensemble 



of 1000 initial conditions with {p n (0)} fulfilling the rela- 
tion in Eq. (|22p with E^ in — 0.1 is used. Note that the 
amount of kinetic energy injected is small compared to 
the energy of the asymmetrical polaron, which amounts 

to E polaron = Y,n\ U n<t>l -2<t>n<t>n-l-{l/1)<t>n] ~ -1-965. 

Notably, after an initial transient where the current den- 
sity grows continually towards larger positive values, it 
reaches an almost stationary regime characterised by 
subdued variations. Most importantly, J stays positive 
throughout the simulation time indicating a positive net 
current density complying with the bias direction induced 
by the asymmetric potential. On the other hand, for re- 
stored reflection symmetry, i.e. U na + n H> U no - n , there 
results zero net current density. 

In summary, we have proposed a transport-rectifying 
mechanism accomplishing directed charge transport in 
the semiclassical Holstein model. We have shown that 
using localised initial conditions, constituted by a po- 
laron solution, together with a nonreversion symmetric 
electron on-site potential suffices for the rectification of 
the current. Crucially, violation of parity prevents the 
existence of pairs of counterpropagating trajectories em- 
anating from an ensemble of unbiased localised initial 
conditions despite the time-reversibility of the underly- 
ing equations of motion. Our numerical simulations of 
the semiclassical Holstein system with an asymmetric 
on-site potential reveal a regime of coherent long-range 
charge transport. Our method can directly be applied 
to induce directed transport in other lattice systems sup- 
porting (mobile) localised solutions. 
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